N = [2:65, 70:5:200]';

sigma = 0.1;

v = asymptoticvbar(sigma, N);

plot(N,v)

omega = 4*pi;
epsilon = sqrt(omega*sigma./(N*pi));


ylim([0 max(v)+0.2]);

hold on;

f = @(x, xdata)fitFunction(x,xdata,sigma);

x = lsqcurvefit(f,1,N,v);
%x = lsqcurvefit(@fitFunction,1,N,v);

plot(N,fitFunction(x,N,sigma),'g');
hline = refline([0 1/15]);
set(hline,'Color','r')
legend('Asymptotic','Homogeneous');
xlabel('N');


hold off;